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I. INTRODUCTION 



We consider phenomenological time reversal invariant violating (TRIV) and parity con- 
serving (PC) interaction, to which we refer as TVPC interaction in order to distinguish it 
from TVPV interactions, which violate both time reversal invariance and parity. TVPC 
interaction was introduced for the first time in paper [l] as a possible explanation of CP- 
violation in ii'^-meson decay. According to CPT theorem, the violation of CP invariance 
implies the TRIV. In spite of the fact that almost all known possible mechanisms of CP 
violation also violate parity, TVPC interactions have been a subject of experimental and 



theoretical studies for decades (see, for example, |2|-|l8| and references therein), because 
they directly manifest phenomena beyond the Standard model. The most experimental con- 
straints for this interactions were obtained by using low energy nuclear physics processes, 
which cover a large variety of nuclear reactions and nuclear decays. There are a number of 
advantages of the search for TRIV in that processes: for example, the possibility of enhance- 
ment of T- violating observables in neutron indused reactions by many orders of a magnitude 



due to complex nuclear structure (see, i.e. paper 10| and references therein), similar to the 



enhancement observed for parity violating effects. Another advantage to be mentioned is the 



19|-|2l|. Then, 



existence of observables which cannot be imitated by final state interactions 
the measurement of non-zero value for these observables directly indicate TRIV, similar to 
the case of neutron electric dipole measurements. 

A promising process for a search for TRIV in nuclear reactions is a measurement of TVPC 

nn 

effects in transmission of polarized neutron through polarized target p|, |4| . These effects can 
be enhanced fl Q by a factor of 10^ and therefore, could be measured at new spallation 
neutron facilities, such as the SNS at the Oak Ridge National Laboratory or the J-SNS at 
J-PARC, Japan. 

However, despite the advantage of the enhancement, complexity of nuclear system makes 
it difficult to directly relate observation of TRIV effects to nucleon TVPC coupling constants. 
Therefore, it is interesting to compare the calculations of TVPC effects in complex nuclei 
with the calculations of these effects in simplest few body systems, which could be useful 
for clarification of infiuence of nuclear structure on values of TVPC effects. Thus, as a 
first step to many body nuclear effects, we study TRIV and parity violating effects in 
one of the simplest available nuclear process, namely elastic neutron-deuteron scattering. 



The calculations of these effects for a specific type of TVPC interaction in a short-range 
approximation show strong dependence of TVPC observables on neutron energy, which 
gives the opportunity to improve existing constrains on TVPC interactions using simple few- 
body system. Therefore, it is desirable to calculate these effects for general case of TVPC 
interactions to clarify this opportunity. 

In this paper we treat TVPC nucleon-nucleon interactions as a perturbation, while non- 
perturbed three-body wave functions are obtained by solving Faddeev equations for realistic 
strong interaction Hamiltonian, based on AV18+U1X interaction model. For description 
of TRIV potentials, we use both meson exchange model and effective field theory (EFT) 
approach. For different cases of symmetry violation we use following notations: "fP for 
TVPC, tf for TVPV, and f for parity violating (PV) cases. 



II. OBSERVABLES 



We consider TVPC effects related to cr„ ■ [p x I]{p ■ I) correlation, where cr„ is the 
neutron spin, / is the target spin and p is the neutron momentum, which can be observed in 
the transmission of polarized neutrons through a aligned (or tensor polarized) target. This 



correlation leads to the difference 



J] between the total neutron cross sections for neutron 



polarized perpendicular to the neutron momentum and averaged over elongated target states 



and neutron spin rotation angle 



Aa^^ = ^Im(/+-/_), (1) 

a, |9j around the axis [p x /](p • J), 

2txN 



MU-f-)- (2) 



dz p 

Here, /+^_ are the zero angle scattering amplitudes for neutrons polarized parallel and anti- 
parallel to the [p X I]{p ■ I) axis respectively, z is the target length, and is a number of 
target nuclei per unit volume. It should be noted that for a non-zero value of this five-fold 
correlation, the spin I must be larger or equal to one, i.e. these TVPC effects require a 
tensor polarized target. Therefore, this correlation cannot be observed in nucleon-nucleon 
scattering. 

The scattering amplitudes can be represented in terms of matrix R which is related to 
scattering matrix S" as ^ = 1 — S*. In partial wave basis, we define Ri>si is = {l'S'\R^\lS) , 



where unprimed and primed parameters correspond to initial and final states, / is an orbital 
angular momentum between neutron and deuteron, 5 is a sum of neutron spin and deuteron 
total angular momentum, and J is the total angular momentum of the neutron-deuteron 
system. Since we are interested in low energy neutrons, one can consider only s, p, and d 
partial waves with mixing only between s and d, and p and p waves. Then, one can write 
the TVPC parameters as 

1 dS^P 1 TT 



N dz 
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-Re 



U„,Z„ ^91^9 ^O'^O ^9'^9 



(3) 
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(4) 



The symmetry violating R -matrix elements can be calculated with a high level of accu- 
racy in Distorted Wave Born Approximation (DWBA) as 



(5) 



where is a neutron-deuteron reduced mass, V^^ is TVPC nucleon-nucleon potential, and 
l^', {l'S')JJ^)^^^ are solutions of 3-body Faddeev equations in configuration space for strong 
interaction Hamiltonian satisfying outgoing(incoming) boundary condition. The factor 
in this expression is introduced to match the i?-matrix definition in the modified spherical 



harmonics convention 



22j with the wave functions expressed in spherical harmonics con- 



vention. The matrix elements of TVPC potential in spherical harmonics convention and 
i?-matrix in modified spherical harmonics convention are antisymmetric under the exchange 
between initial and final states. 

For calculations of wave functions, we used jj-coupling scheme instead of IS coupling 
scheme. We can relate i?-matrix elements in IS coupling scheme to jj-coupling scheme using 
unitary transformation (see, for example 23l |) 



[Sk ® Ja 



)s\jJ^ = ^\[jx ® {ly ® Sfc)jJjjJ 



3y 



ly Sk iy 



x(-i)>+i.-^(_l)'.+«^+>+^[(2jj^ + l)(25+ 1)]^ <( ' )>,(6) 

jx J S 



where jx is a spin (total angular momentum of the target, jx = 1 for the deutron), Sk is a 
spin of the projectile [sk = \ for the neutron). In IS coupling scheme spins of the projectile 
and the target are added giving partial spin S to which relative projectile-target angular 
momentum (/) is added to obtain total angular momentum (J) of the system. In contrary, 
in jj-coupling scheme, the relative projectile-target angular momentum {I) is added to the 
projectile spin (s^) giving intermediate angular momentum (j) before coupling it with target 
spin {jx) in order to obtain total angular momentum (J) of the system. 



III. TIME REVERSAL VIOLATING PARITY CONSERVING POTENTIALS 



The most general form of time reversal violating and parity conserving part of nucleon- 
nucleon Hamiltonian in the first order of relative nucleon momentum can be written as 



P 

rriN 

P 



+ {9b{r) + fi'6(r)ri ■ T2 + 97{r)T^2 + fi'8('^)^+) en ' (^2r 

+ {99{r) + 9io{r)Ti ■ T2 + 9u{r)T^2 + 9i2{r)T+) 

^ ( - P , - P 2„ p 
X r ■ CTi ■ 0-2 + r ■ (72 ■ (Ti- -r ■ cti ■ (T2 

+ i9i3ir) + 9u{r)Ti ■ T2 + 9ih{r)T^2 + 9i&{.'r)r+) 

f ^ ^ ^ P ^ ,^ P P P 
X r ■ (TiT ■ (T2r ■ -(r ■ cti ■ <T2 + r ■ cri ■ (T2 + r ■ (T2 ■ (Ti 

\ rriN 5 rriN tun rnN 

+^i7(r)r_f ■ (cr^ x ^) + ^isWr^f • x ^), (7) 
rriN- rriN 

where exact form of 9i{r) depends on the details of a particular theory of TVPC 

One should note, that pions, being spin zero particles, do not contribute to TVPC on- 



shell interaction 25|. Therefore to describe TVPC nucleon-nucleon interactions in meson 
exchange potential model, by assuming CPT conservation, one should consider contribution 
from heavier mesons: p(770), /^( J^^) = and /ii(1170), /^( J^^) = 0-(l+-) (see, 

for example . Il2| and references therein). For example, Lagrangians for p and hi are 



= -9pN{^,p^'^ - ^a.^d-^p'^nr^'N - 9hNi'i,KN, (8) 
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where we neglected terms, such as N'^^d^h^N ^ which are small at low energy, and g and 
g represent strong and TVPC meson nucleon couplings respectively. Then, one can obtain 
TVPC potentials 



l\e~x 

\ME 



.rfp 9h9hmf^^ r w P . , P /,^^ 

= — ^ Yi{mhr){ai- -—CT2-r + cr2- -^(Ti-r), (10) 

where Yi{x) = (1 + ^)^-^, Xa = ruar. Comparing these potentials with eq. ([7]), one can see 



that in this model, all gi{r) = 0, except for 



2gh9h 

rriN 



where we introduced dimensional constants and scalar function /^'^(/i) = |^Y'i(/ir), so 
that gn{r) can be written as 

9n{r) = J2cZfta{r). (12) 

a 

If we include iso-vector J = 1 ai and 6i mesons, which masses are close to the value of hi 
mass, functions g^ and giQ will also contribute to TVPC potential. 

In EFT approach we consider eq.© as a leading order of TVPC potential. Final result 
must not depend on the particular form of the gn{r) functions as long as they are localized, 
like delta function or its derivative. Therefore, using EFT one can estimate the contribution 
of each term of the potential substituting (7„(r) functions by the corresponding products of 
low energy constants (LECs) and Yukawa functions Yi{fir). The mass scale fi represents a 
regularization scale of EFT. For example, at low energy we can assume that /j, ~ for 
pionless EFT approach. 



IV. RESULTS AND DISCUSSIONS 



For calculation of TRIV amplitudes in DWBA approach, we used the non-perturbed 
(time reversal invariance conserving) 3-body wave functions for neutron-deuteron scattering 
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obtained by solving Faddeev equations (also often called Kowalski-Noyes equations) in con- 
figuration space 0, Is?]. The detailed procedure for these calculations is described in our 
papers 23|, |28|. As previously we employed Ayi8 nucleon-nucleon potential in conjunction 
with UIX three-nucleon force. Obtained contribution of each TVPC operator from eq. ([7]) 
to the matrix element of eq.(|5]) is summarized in Tables [I] and [ITl representing real and 
imaginary parts respectively. The matrix elements are evaluated using jj coupling scheme 
for the neutron-deutron center of mass energy equal Ecm = 100 keV and the regulariza- 
tion scale set by /i = rup, which is equal to the mass of the lightest meson contributing to 
TVPC interaction. It should be noted that each matrix element presented in these tables 
contains a sum of contributions from different Faddeev components of wave functions with 
a large number of partial waves. Therefore, the values of the matrix elements are strongly 
dependent on detailed behavior of exact wave functions. However, in spite of the fact of the 
possibility of a numerical suppression of matrix elements for some operators, the calculated 
values are stable enough to be used for estimations of TRIV effects. 

The contributions from each TVPC operator to the difference of scattering amplitudes 
/+^_ are summarized in Table IIIH where we distinguish three columns representing result 
with a different choice of the characteristic mass scale for gi{r) functions. Thus, for example, 
the column ^y- corresponds the description of TVPC potential in pionless EFT. (As it was 
mentioned above, vr-meson exchange cannot lead to TVPC interaction.) 

It should be noted that in spite of the fact that all results of calculations are presented 
only for neutron energy Ecm = 100 keV, they can be easily extrapolated for any value of 
neutron energy below 1 MeV since they have a simple dependence on neutron energy ii^ as: 

Re— ^ Ve, Im— E. (13) 

p p 

This is because these TVPC observables are the result of a mixing of initial and final p- waves, 
or s- and ci- waves, by TVPC interactions in scattering amplitudes. 

To have insights into the structure of TVPC scattering amplitudes, one can compare 
them with strong, PV, and TVPV amplitudes at the same energy E^m = 100 keV {p = 
0.567 X 10-^ fm"^), and /i = m^, which corresponds to pionless EFT potential. Then, 
strong s-wave scattering amplitude f^^ is 

£st 

y = (72.2 + i24.8) /m^ (14) 



TABLE I. Representative contribution of each TVPC potential term to the real part of the matrix 
element (Re ''''^^ ^ij — ^^^"''^''^^ ). Results are presented using jj-coupling scheme for wave functions 
obtained using ^41/18 + UIX interaction at Ecm = 100 keV. For all operators a scalar function 
-^Yi{mpr) has been used. Operators 3,7,11 and 15 are null due to the isospin selection rules. All 
data are in fm?. 
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giving the total cross section atot = = 3.11 b. PV difference of scattering 

amplitudes in EFT is 23 1 

1 A/^(/i = m^) 



[(-1.93 ■ ■ • 2.42) + i(-0.22 ■ ■ • 0.67)] fm^ 

rriNCn P 

The difference of TRIV amplitudes with parity violation in EFT is 
1 A/^^(/i = m.) 

and for TVPC ones is 

1 A/7^(/i = m.) 



(15) 



[(-1.63 ■ ■ ■ 0.66) + i(-0.063 ■ ■ ■ 0.22)] fm\ (16) 



tunC^ 



TP 



P 



[(-0.01 ■ ■ ■ 0.03) + i(-0.0013 ■ ■ ■ 0.0004)] /ml (17) 
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TABLE II. Same as in Table U but for the imaginary part of the matrix element 
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Here, C^, C^^, and C^^ are low energy constants for PV, TVPV, and TVPC interactions, 
correspondingly, with a separated factor I/itln to match dimensions of the final result and 
retain dimension of LECs C„ in [/m]. The range of values for real and imaginary parts 
of PV and TRIV amplitudes is defined by a value of possible contribution from each PV, 
TVPV or TVPC operator. 

Aforementioned amplitudes follow the simple kinematic rule of the suppression ~ {pRnuc) 
for an additional p-wave involved in PV and TVPV amplitudes, and ~ {pRnuc)'^ for two p- 
waves or one (i-wave for the case of TVPC amplitudes. (Here, Rnuc is an effective range of 
strong interaction, which leads to ~ (pRnuc) ~ 0.1 for neutron energy Ecm = 100 keV.) In 
addition to this kinematic factor, TVPC scattering amplitudes are suppressed, as compared 
to PV or TVPV ones, by a factor — ~ 0.1, which results from an extra momentum 
dependence of all operators in TVPC potential. By increasing neutron energy, one can 
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TABLE III. Difference of scattering amplitudes, (/-f — f-)/p for TVPC potential from each opera- 
tors and mass scales at Ecm = 100 keV. Note that pion mass scale does not corresponds to physical 
meson exchange potential. Operators with n = 3, 7, 11, 15 are equal to zero due to isospin selection 
rules. All data are in frm?. 
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easily increase the kinematic factor up to one. Then, the only suppression of TVPC matrix 
elements in the amplitude will be left due to ~ 0.1 . It should be noted that this 

^ ^ . . . . mjv 



suppression factor is well known |8|-|10 



18| for TVPC matrix elements in nuclei. 
It is noteworthy that our calculations are in good agreement with results Q, obtained 
using zero range force approximation for calculations of TVPC effects m n — d scattering. 
For example, using eq.(8) of paper one can obtain for Ecm = 100 keV 

-J— = ^/(0.0004 + zO.0013) fm^, (18) 

where g/ is unknown TVPC nucleon-nucleon coupling constant. 

The results of table UTTl can also be used to express TVPC parameters in terms of meson 
exchange model. Since TVPC meson exchange model does not allow pion exchanges, the 

10 



lightest mesons to be considered are p and hi mesons. Then, assuming only contributions 
from these mesons, one can obtain for Ecm — 100 keV 



Finally, we conclude that neutron-dcuteron scattering is a promising process to improve 
current experimental constraints on TVPC interactions. The TVPC observables can be 
large enough to be measured at neutron energy of hundreds of keV due to strong energy 
dependence. On the other hand they can be precisely calculated, providing the possibility 
to extract the TVPC nucleon coupling constants from the experiment. 
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